Three dimensional simulation of fluid flow in X‐ray CT images of porous media     

A. Al‐Omari  E. Masad 《国际地质力学数值与分析法杂志》2004,28(13):1327-1360

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

A comparison of theoretical Sv emission line strengths with extreme ultraviolet observations of a sunspot     

F. P. Keenan  J. G. Doyle 《Solar physics》1990,128(2):345-351

Electron impact excitation rates for transitions in the S v ion, calculated with theR-matrix code, are used to derive the electron temperature sensitive emission line ratiosR ₁ =I(854.8 Å)/I(786.9 Å),R ₂ =I(852.2 Å)/I(786.9 Å),R ₃ =I(849.2 Å)/I(786.9 Å), andR ₄ =I(1199.1 Å)/I(786.9 Å), which are found to be significantly different from previous estimates. A comparison of the present results with observational data for a sunspot obtained with the Harvard S-055 spectrometer on boardSkylab reveals generally good agreement between theory and experiment, except in the case ofR ₁, which is probably due to blending in the 854.8 Å feature. The possible effects of Lyman continuum absorption on the observed line ratios is briefly discussed.  相似文献   

The dust coma of periodic comet Churyumov-Gerasimenko (1982 VIII)     

M.S. Hanner  E. Tedesco  A.T. Tokunaga  G.J. Veeder  D.F. Lester  F.C. Witteborn  J.D. Bregman  J. Gradie  L. Lebofsky 《Icarus》1985,64(1):11-19

The dust coma of Comet P/Churyumov-Gerasimenko was monitored in the infrared (1–20 μm) from September 1982 to March 1983. Maximum dust production rate of ～2 × 10⁵ g/sec occured in December, 1 month postperihelion. The ratio of dust/gas production was higher than that in other short-period comets. No silicate feature was visible in the 8- to 13-μm spectrum on 23 October. The mean geometric albedo of the grains was ～0.04 at 1.25 μm and ～0.05 at 2.2 μm.  相似文献   

The geochemistry of lead in rivers,estuaries and the continental shelf of the Southeastern United States     

Herbert L Windom  Ralph G Smith  Masaru Maeda 《Marine Chemistry》1985,17(1):43-56

Transport by southeastern rivers has insignificant influence on the lead concentration of southeastern shelf waters. If the rate of mobilization and transport of lead by these rivers is representative of uncontaminated fluvial transport during the Pleistocene, only about 5% of the prehistoric output of lead to the North Atlantic can be accounted for by river input.Lead concentrations in southeastern shelf waters are similar to those of North Atlantic Water from the upper 1000 m, which is probably similar to the intrusion source water for the shelf. Atmospheric inputs to the shelf of the same magnitude as observed for the Western North Atlantic are difficult to reconcile given the residence time of shelf waters and their lead concentration unless the rate of loss of lead to shelf sediments is about the same as the atmospheric flux.  相似文献   

Consideration of time-correlated errors in a Kalman filter applicable to GNSS   总被引：1，自引：0，他引：1  

M. G. Petovello  K. O’Keefe  G. Lachapelle  M. E. Cannon 《Journal of Geodesy》2009,83(1):51-56

An algorithm for considering time-correlated errors in a Kalman filter is presented. The algorithm differs from previous implementations in that it does not suffer from numerical problems; does not contain inherent time latency or require reinterpretation of Kalman filter parameters, and gives full consideration to additive white noise that is often still present but ignored in previous implementations. Simulation results indicate that the application of the new algorithm yields more realistic and therefore useful state and covariance information than the standard implementation. Results from a field test of the algorithm applied to the problem of kinematic differential GPS demonstrate that the algorithm provides slightly pessimistic covariance estimates whereas the standard Kalman filter provides optimistic covariance estimates.  相似文献   

Pressure-induced incipient amorphization of α-quartz and transition to coesite in an eclogite from Antarctica: a first record and some consequences     

R. PALMERI  M. L. FREZZOTTI  G. GODARD  R. J. DAVIES 《Journal of Metamorphic Geology》2009,27(9):685-705

Monocrystalline quartz inclusions in garnet and omphacite from various eclogite samples from the Lanterman Range (Northern Victoria Land, Antarctica) have been investigated by cathodoluminescence (CL), Raman spectroscopy and imaging, and in situ X‐ray (XR) microdiffraction using the synchrotron. A few inclusions, with a clear‐to‐opalescent lustre, show ‘anomalous’ Raman spectra characterized by weak α‐quartz modes, the broadening of the main α‐quartz peak at 465 cm^?1, and additional vibrations at 480–485, 520–523 and 608 cm^?1. CL and Raman imaging indicate that this ‘anomalous’α‐quartz occurs as relicts within ordinary α‐quartz, and that it was preserved in the internal parts of small quartz inclusions. XR diffraction circular patterns display irregular and broad α‐quartz spots, some of which show an anomalous d‐spacing tightening of ～2%. They also show some very weak, hazy clouds that have d‐spacing compatible with coesite but not with α‐quartz. Raman spectrometry and XR microdiffraction characterize the anomalies with respect to α‐quartz as (i) a pressure‐induced disordering and incipient amorphization, mainly revealed by the 480–485 and 608‐cm^?1 Raman bands, together with (ii) a lattice densification, evidenced by d‐spacing tightening; (iii) the cryptic development of coesite, 520–523 cm^?1 being the main Raman peak of coesite and (iv) Brazil micro‐twinning. This ‘anomalous’α‐quartz represents the first example of pressure‐induced incipient amorphization of a metastable phase in a crustal rock. This issue is really surprising because pressure‐induced amorphization of metastable α‐quartz, observed in impactites and known to occur between 15 and 32 GPa during ultrahigh‐pressure (UHP) experiments at room temperature, is in principle irrelevant under normal geological P–T conditions. A shock (due to a seism?) or a local overpressure at the inclusion scale (due to expansion mismatch between quartz and its host mineral) seem the only geological mechanisms that can produce such incipient amorphization in crustal rocks. This discovery throws new light on the modality of the quartz‐coesite transition and on the pressure regimes (non‐lithostatic v. lithostatic) during high‐pressure/UHP metamorphism. In particular, incipient amorphization of quartz could favour the quartz‐coesite transition, or allow the growth of metastable coesite, as already experimentally observed.  相似文献   

Interplay between erosion and tectonics in the Western Alps     

Marco G. Malus  Giovanni Vezzoli 《地学学报》2006,18(2):104-108

Bedrock fission‐track analysis, high‐resolution petrography and heavy mineral analyses of sediments are used to investigate the relationships between erosion and tectonics in the Western Alps. Along the Aosta Valley cross‐section, exhumation rates based on fission‐track data are higher in the fault‐bounded western block than in the eastern block (0.4–1.5 vs. 0.1–0.3 mm yr⁻¹). Erosion rates based on the analysis of bed‐load in the Dora Baltea drainage display the same pattern and have similar magnitudes in the relative sub‐basins (0.4–0.7 vs. 0.04–0.08 mm yr⁻¹). Results highlight that climate, relief and lithology are not the controlling factors of erosion in the Western Alps. The main driving force behind erosion is instead tectonics that causes the differential upward motion of crustal blocks.  相似文献   

Observational constraints on the radio and γ-ray emission regions of PSR B1055−52     

H. G. Wang  G. J. Qiao  R. X. Xu  Y. Liu 《Monthly notices of the Royal Astronomical Society》2006,366(3):945-952

  相似文献   

Cave ice of Lake Baikal as an indicator of climatic changes     

E. V. Trofimova 《Doklady Earth Sciences》2006,410(1):1087-1090

  相似文献   

P‐wave transmission across fractures with nonlinear deformational behaviour     

X. B. Zhao  J. Zhao  J. G. Cai 《国际地质力学数值与分析法杂志》2006,30(11):1097-1112

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献